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We further explore the consequences of treating the X(3872) meson as a tetraquark bound state by 
analyzing its one-photon decay X 7 + J/xp in the framework of our approach developed in previous 
papers which incorporates quark confinement in an effective way. To introduce electromagnetism we 
gauge a nonlocal effective Lagrangian describing the interaction of the X(3872) meson with its four 
constituent quarks by using the P-exponential path-independent formalism. We calculate the matrix 
element of the transition X — >■ 7 -I- J/i/) and prove its gauge invariance. We evaluate the X ^> 7 -f J/i/j 
decay width and the longitudinal/transverse composition of the J/i/) in this decay. For a reasonable 
value of the size parameter of the X(3872) meson we find consistency with the available experimental 
data. We also calculate the helicity and multipole amplitudes of the process, and describe how they 
can be obtained from the covariant transition amplitude by covariant projection. 
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I. INTRODUCTION 

This paper is a direct continuation of our previous work ll| where we have analyzed the strong decays of the 
charmonium-like state X(3872) in the framework of our relativistic constituent quark model which includes infrared 
confinement in an effective way [2|. In our approach the X(3872) meson is interpreted as a tetraquark state with the 
quantum numbers J^*^ — 1^^ as in In this paper we analyze the one-photon decay X -> 7 + J/0 in the same 
tetraquark picture. The electromagnetic interaction is incorporated into our relativistic nonlocal effective Lagrangian 
in a gauge invariant way using the P-exponential path-independent formalism. 

We begin by collecting the experimental data relevant for our purposes. A narrow charmonium-like state X(3872) 
was observed in 2003 in the exclusive decay process B'^ — >■ i^='=7r+7r~ J/0 0]. The X(3872) decays into ■n^'n:~ J/il) and 
has a mass of mx = 3872.0 ± 0.6(stat) ± 0.5(syst) very close to the m^jo -I- m£i*a = 3871.81 ± 0.25 mass threshold 
0. Its width was found to be less than 2.3 MeV at 90% confidence level. The state was confirmed in B-decays by 
the BABAR experiment [gJ and in pp production by the Tevatron experiments CDF Q and D0 Q . The most precise 
measurement up to now was done in [0| with mx = 3871.61 ± 0.16 ± 0.19. The new average mass given in Q is 

mx ^ 3871.51 ± 0.22 MeV. (1) 
The Belle Collaboration has reported fTol evidence for the decay modes X(3872) — 7-f J/0 and to AT — > tt+tt^tt^ J/0: 



B{B -> XK) ■ B{X ^ 7 + J/iP) = (1.8 ± 0.6 (stat) ± 0.1 (syst)) x 10"^ 
r(A ^ 7 + J/yj) 



B{X 7r+7r-7rOj/0) 

B{X TT+TT-J/lP) 



= 0.14 ±0.05, 

= 1.0 ±0.4 (stat) ±0.3 (syst). (2) 



These observations imply strong isospin violation because the three-pion decay proceeds via an intermediate uj 
meson with isospin whereas the two-pion decay proceeds via the intermediate p meson with isospin 1. It is evident 
that the two-pion decay via the intermediate p meson is very difficult to explain by using an interpretation of the 
A (3872) as a simple cc charmonium state with isospin 0. 

In an analysis of 5+ J/%1} 7 decays, the BABAR Collaboration [ll| found evidence for the radiative decay 
A (3872) — )■ 7 + J/0 with a statistical significance of 3.4(t. They reported the following values for the product of 
branching fractions 

B{B+ XK+) ■ B{X ^ 7 + J/V-) = (3.3 ± 1.0 (stat) ± 0.3 (syst)) x 10"^ . (3) 

The Belle Collaboration reported [12] the first observation of a near-threshold enhancement in the D^D^-k^ system 
from B D"D^n°K. The enhancement peaks at a mass of M = 3875.2 ± 0.7^5 j ± 0.8 MeV. The branching fraction 
for events in the peak is 

B{B ^ D"D\°K) = (1.22 ± 0.3lj:°:23) x lO"" . (4) 
All available experimental data up to 2007 were analyzed in [13j. The authors found that [l^ 

B{B+ ^ AA:+) l.SOtail X lO^"* , 

r(A ^ 7 + j/v-) 



r(A ^ TT+TT-J/tP) 



= 0.22 ±0.06. (5) 



The BABAR Collaboration found evidence for the decays A ^ 7 + J/0 and A — >■ 7 + ■0(25') in their data sample of 
the B — )■ ccjK decays. The measured products of branching fractions are [l3 | 

B{B^ XK^) ■ B{X ^ 7 + J/ij) = (2.8 ± 0.8 (stat) ± 0.1 (syst)) x 10"^ , 
B{B^ XK^) ■ B{X ^ 7 + ^A(25)) = (9.5 ± 2.7 (stat) ± 0.6 (syst)) x 10"^. (6) 
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There have been many theoretical attempts to unravel the structure of the X(3872) and its decays. Many of the 
theoretical predictions for the decay X(3872) — > 7 + J/'tp published up to now are very model dependent. We mention 
some of them in turn. 

All possible ID and 2P cc assignments for the X(3872) were considered in [15]. The authors obtained E\ radiative 
widths for decays into charmonium cc states as well as for some strong decays taking the experimental mass as input. 
The conclusion was that many of the possible J^'-' assignments can be eliminated due to the smallness of the observed 
total width. The suggestion was that radiative transitions could be used to test the remaining 

jPC 

assignments. 

Some tests of the hypothesis that the X(3872) is a weakly bound D^D^* molecule state were suggested in ^6j]. It 
was proposed that measuring the 37r J/0, 7 + J/i/;, 7 + ?/;', KK*, and np decay modes of the X will serve as a definitive 
diagnostic tool to confirm or to rule out the molecule hypothesis. 

Assuming that the X(3872) state has the structure {D°D° * - D° *D°)/V2 with quantum numbers J^"^ = the 
X(3872) — > 7 + J/0 decay width was calculated using a phenomenological Lagrangian approach (T^]. The calculated 
value of the radiative decay width varied from 125 KeV to 250 KeV depending on the model parameters. 

QCD sum rules were used in [l^ to calculate the width of the radiative decay of the meson X(3872), which was 
assumed to be a mixture between charmonium and exotic molecular [cg][gc] states with J^^ = In a small range 
for the values of the mixing angle, one obtains 



+ 0.19 ±0.13. (7) 



T{X ^ J/0 7r^ 



Our paper is organized as follows. In Sec. II we gauge a nonlocal effective Lagrangian describing the interaction of 
the X(3872) meson with its constituent quarks by using the P-exponential path-independent formalism developed in 
[20| . In Sec. Ill we calculate the matrix element of the radiative transition X — >• 7 + J/0 and prove its gauge 
invariance analytically. In Sec. IV we present the results of our numerical analysis. First, we check numerically that 
the final amplitude is gauge invariant. Second, we introduce infrared confinement as was done in our previous papers 
Refs. 1, 2] and evaluate the X — > 7 + J/0 decay width. Finally, in Sec. V we summarize our results. In an Appendix 
we describe how the two helicity or the two multipole amplitudes of the process can be obtained from the gauge 
invariant transition amplitude by covariant projection. 



II. THEORETICAL FRAMEWORK 



The effective interaction Lagrangians describing the coupling of the charmonium-like meson such as the A'(3872) 
to four quarks, and the coupling of the charmonium J/0 state to its two constituent quarks are written in the form 
(see Ref. [l|) 



^int ^ gxXg^ix)- Jf^^ (x) + gj/^ J/0^ (x) • J^^^ {x) {q^u,d). (8) 
The nonlocal interpolating quark currents read 



Jx, (x) = / dxi .. . dxi ^ ( a; - ^ WiXi j $x ^(a^j - Xj Y^ 
"J J \ i=i J i<j 

X -^eabcSdec { fea (a;4)C7^Ch (a;i )] (2:3 C'ce (x2 )] + (7^ o 7^) |, 



nic Wc 

Wi = W2 = —, r = — , ^3 = ^4 



2{mq + rric) 2 2(mq + rric) 2 

The matrix C = 7^72 is related to the charge conjugation matrix: C = Ct = C^^ = -C^, CT^C^^ = ±T, (" + " 
for T = S, P, A and " — " for F — V, T) . We follow Q and take the tetraquark state to be a linear superposition of 
the X„ and Xd states according to 
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Xi = Xiow = Xu COS 6 + Xd sin 6, 
Xh = Xhigh = -Xu sin e + Xd cos 6*. (10) 

The coupling constant gx in Eq. ([S]) will be determined from the compositeness condition Zh — (see e.g. 
Refs. pll. [23]1. The compositeness condition requires that the renormalization constant Zh of the elementary meson 
X is set to zero, i.e. 

ZH = l-Il-'H{VH=nil)^Q, (11) 

where IVxip^) is the scalar part of the meson mass operator and the prime stands for the derivative w.r.t. p^j. For 
the spin one states X(3872) and J/i) the compositeness condition reads 

nv{p')^\{g,.~^-^^K^{p). (12) 

The X meson mass operator can be calculated from the self-energy three-loop sunrise-type diagram with four quark- 
antiquark propagators. The calculation is described in more detail in Ref. [l|. 

As in the case of baryons composed of three quarks it is convenient to transform to Jacobi coordinates in the 
integrals of Eq. In the case of four quarks one has 



^^^"^^ — — — 271"^^ ^ — 2 — ^"^^^^ 



2wi + 1113 + W4 W3 — W4 W3 + W4 \ ^ 

= ^ ^ — + ^^p' = " + ^ 

3 

= — ^ — P2-^^P3 = x+x.^3,ft' 

3 

Wi - W2 Wi + W2 + 2W3 Wi +W2 v-^ „, 

^^^^ 27^^^ 271 2 = ^ + '^"jT'i, (13) 

4 3 

where x = ^ XiWi and ^ (2:^ — Xj)^ = ^ . The inverse transformation reads 

i=l l<*<i<4 j=l 

Pi = V2{xi - X2), P2 ^ V2ix3 - P3 = Xi + X2 - X3 - X4. 

In the case of two quarks as e.g. in the J/0 case one has 

11 

yi = y + -^p, y2 = y- -^p- (14) 

One then has 

= j dp<i>x{p'^)J4q{xu---,X4), 

J^g{xi,...,X4) = -^eabcSdec^ fea (a;4)C7^C,, (xi )] [^d (^3 )7''C'Ce (x2 )] + (7^ ^ 7'') |, 
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Jj/^iy) = ]dp<^j/^{p^)J!^,{yuy2), J^,{yi,y2) = ca{yih''ca{y2), (15) 

where dp — dpidp2dp3 and p^ — Pi + P2 ^ Pa- The Jacobian is absorbed into the coupling gx- 

The gauge invariant interaction of a bound quark state with the electromagnetic field has been described in some 
detail in Ref. [l9| . For comprehensive purposes we recall some of the key points of the gauging process. Since the 
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X{3872) and J/0 mesons are neutral mesons we will discuss the charged quarks only. The free Lagrangian of quarks 
is gauged in the standard manner by using minimal substitution: 

d^q-^ {d" -iegA'')q, d^q ^ {d" + iegAi')q, (16) 

where Cg is the quark's charge (en = | e, — ^'^'^^ etc.). Minimal substitution gives us the first piece of the 
electromagnetic interaction Lagrangian 

'^(^) ^E^^^mW-/,"!^). J!:ix)^qix)rqix). (17) 

Q 

In order to guarantee gauge invariance of the nonlocal strong interaction Lagrangian, one multiplies each quark 
field q{xi) in the relevant quark current J^{x) given by Eq. (jl5p by a gauge field exponential according to 

q{x,) ^ e-'^-'^^-^^P^ q{x,), q{x{) ^ e^^'^^^-^'^^ q(x,), 

Xi 

I{x,,x,P) - / dz^A^(z). (18) 



where P is the path taken from a; to x^. It is readily seen that the full Lagrangian Eq. ([8]) is invariant under the local 
gauge transformations 

q{x,) ^ e'^'/(-')g(x,), q{x,) ^ e"- '/(-')g(:r,), 

A^{z)^ A^'{z) + ^^'f{z), so that I{x,,x,P) ^ I{x,,x,P) + f{x,) - f{x). (19) 

The second term of the electromagnetic interaction Lagrangian >Cf,™.2 arises when one expands the gauge exponential 
in powers of up to the order of perturbation theory that one is considering. Superficially the results appear to 
depend on the path P which connects the endpoints in the path integral in Eq (jlSp . However, one needs to know 
only derivatives of the path integrals when doing the perturbative expansion. One can make use of the formalism 
developed in [13] which is based on the path-independent definition of the derivative of I{x, y, P): 

Jim dx^'j^I{x,y,P) ^ \im [I{x + dx,y,P')~Iix,y,P)], (20) 

where the path P' is obtained from P by shifting the endpoint x by dx. Use of the definition ((20|) leads to the key 
rule 



d 



Iix,y,P) = A^{x) (21) 



which states that the derivative of the path integral I{x, y, P) does not depend on the path P originally used in 
the definition. The nonminimal substitution (jl8p is therefore completely equivalent to the minimal prescription as is 
evident from the identities (I20p or (j2ip . The method of deri ving Feynman rules for the nonlocal coupling of hadrons 
to photons and quarks was worked out before in Refs. [l^ |2(| and will be discussed in the next section where we 
apply the formalism to the physical processes considered in this paper. 
Expanding the Lagrangian up to the first order in A^^ one obtains 

^'int''^\^) = 9x Xgf,{x) ■ J^^_^^{x) + gj/^ JAAp(x) • Jj^^_^Jx) (q = u,d), 
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In order to use the key rule Eq. (j2ip we take the Fourier-transforms for the vertex functions $ and quark fields q 
One then writes down 



(27r)4 ^'^-^ J (27r 



-i— 1 



i?;^ = 9; + ^p^ (24) 



pi 

where 



Wl = CiiPi - C21P2 - C31P3 + C41P4, 
W2 = C12P1 - C22P2 - C32P3 + C42P4, 

W3 = C13P1 - C23P2 " C33P3 + C43P4, 
P=l(pi+P2). (25) 
Finally, we employ a convenient identity which was proven in [lol ] . The identity reads 

1 

FiD%)IT = 1 dr^^'(ri?^^ - (1 - r)a.f ) c,, (a;^.A.(:.,) - 2 z c.JA.(x,)) + Fi-oj^)!^ . (26) 


The identity holds for any function F{z) that is analytical at 2; = 0. 
One obtains 



Jx,-cmix) = n/ J d^y Jiqixi, . . . ,X4) Ap{y) E'^{x;xi, . . . ,X4,y), 



(27) 



4 = 1 

4 . ,4 



i;^(x;xi,...,x4,y) = J 



^4 _ 

^-ipi(x-xi)+ip2{x~X2)+ip3{x-X3)-ip4{x-X4)-ir{x-y) j^p^^p^ 



(2^)4^ 



'^i = (cjj-r + 2wj), (i = 1,. . . ,4; j = 1, . . .,3), 

2^11 = T [ciir + wi)^ + (1 - r) + ^2 + '^s: 

Zi2 = {ctir + wi)^ + r (ci2r + a;2)^ + (1 - r) + ^3, 

= (cji?" + LOif + {ci2r + + T {ctsr + uj^f + (1 - r) loI- 
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(2/) 



d'^yi I d'^V2 / d^zJ2q(yi,y2)Ap(z)E'lj^^{y;yi,y2,z), 



(28) 



Ej,Jy;yi,y2,z) = 



(27r)4 J (27r)4 J (27r) 



? -»Pi(yi-y)+«P2fe2-y)+«9(^-y) EPj,^{pi,p2,q) , 



"j/^ipi,P2,q) = ej dT[- <i>'j/4~z^) E - <i>'j/4-z+) , 



^ t{pT ^q) - - t)p'^ , l^=pzplq^ p^\{pi+P2). 

For calculational convenience we will choose a simple Gaussian form for the vertex function ^^). The minus 

sign in the argument of the Gaussian function is chosen to emphasize that we are working in Minkowski space. One 
has 



$x(-f^') =cxp(l]VA|) 



(29) 



where the parameter Kx characterizes the size of the X meson. Since fi^ turns into — in Euclidean space the form 
(l29l) has the appropriate fall-off behavior in the Euclidean region. We emphasize that any choice for $x is appropriate 
as long as it falls off sufhciently fast in the ultraviolet region of Euclidean space to render the corresponding Feynman 
diagrams ultraviolet finite. As mentioned before we shall choose a Gaussian form for <^x in our numerical calculation 
for calculational convenience. 



III. MATRIX ELEMENT FOR THE DECAY X ^ f + J/tp 

The matrix element of the decay A" (3872) — > ^ + J/ij] can be calculated from the Feynman diagrams shown in Fig.[T] 
The invariant matrix element for the decay is given by 





(c) (d) 

FIG. 1. Feynman diagrams describing the decay A — >■ 7 + J/tp. 



MiXgip) ^ J/V'(gi)7('Z2)) = *(2^)45(4)(p_qi -52)4 ePs''j/^T^p,{qi,q2) 



(30) 
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where 



T^pu (gi , 92 ) = ^MP'' ' 92 ) , 

i—a,b.c,d 



X itr 



75'S'c(fci)7i"S'c(fci + qiht,Sgik2)jpSq{k2 + 92) - (75 ^ 7^1) 



K. 



{ki + \ qif + \{k2 + \ q2f + \ (wqqi - Wcq2) 



2 _ 1 

a — 2 V'^l ^ 2 yi'' T 2 T 2 y2y T 4 v'"9'J 



^ipl = 6 72 5x 5.W y (-27r^4- ^ ( ^ (^2 + \qif) Expijpi,---, Pi, r) 



X jtr 



{2t:)H J {2ny 
lbSq{ki)"ff^Sc{k2)^uSc{k2 + qi) - (75 ^ 7m) 



Pi = ^2, P2 = k2 + qi, Ps^Pi^-ki, r = -52 , 



X itr 



gxgj,^e.jj^Jj^<^x[-K^)'^j/^[-{k2 + q2+2 
l5Sq{ki)"ff^Sc{k2)^pSc{k2 + q2)li^Sc{k2 +p) ~ (75 ^ 7/^) 



91 J 



d^ki f d^k2 



X 



X itr 



Pi 



(27r)4i J (27r)4i 
liiSq{ki)-f^Sc{k2)luSc{k2 +p) - (75 7p) 
-k2~p, P2 = -^2, q = -q2- 



Kc)Ej/^p{pi,P2,q) 



We have analytically checked on the gauge invariance of the unintegrated transition matrix element by contraction 
with the photon momentum q2 which yields qt^Tf^p^iqi, 92) = using the identities 

S{k2) i2 S{k2 + 92) = S{k2 + 92) - S{k2) , 



dr $'(-r a - (1 - r) 6) (a - 6) = $(-6) - $(-a). 



IV. NUMERICAL RESULTS 



The evaluation of the loop integrals in Eq. ([50)1 proceeds as described in our previous paper [l| . If one takes the 
on-mass shell conditions into account 



£xPp = 0, e^j/^qi^ = 0, ePq2p = 
one can write down five seemingly independent Lorentz structures 



(31) 



Tppu{qi,q2) = £q2fiiyp{qi ' 92) Wl + Sq^q^^pqip W2 + eq^q^ppq2u W3 + Eq^q^p^qip W4 + Eq^pypiqi ■ ^2) W5 . (32) 

Using the gauge invariance condition 

q2T^.p. = [qi ■ q2)eq,q,pAWi + W^s) - (33) 
one has W4 — —W5 which reduces the set of independent covariants to four: 

Tpp„{qi,q2) = (<?1 • 92) £q2f^l^p ^1 + £gig2l^p<?lM ^2 + £qiq2p.pq2l^ W3 + (Sq^q^p^qip - {qi ■ (?2)£gipi.p) W4 . (34) 



9 



The gauge invariance condition W4 = —W5 provides for a numerical check on the gauge invariance of our calculation 
as described further on. 

However, there are two nontrivial relations among the four covariants which can be derived by noting [23l | that the 
tensor 

K f 41^5] = 9p-i^i£iy2i^3'^4'^5 + cycl.(j/ii/2J^3i^4i^5) (35) 

vanishes in four dimensions since it is totally antisymmetric in the five indices (1^1,1/2 ,1^3,1^ a, 1^5)- Upon contraction 
with 91(71^92^ and 92 91^92^ finds (between polarization vectors) 

(ll£q2tiup + £gig2fp'?lM ^" (^^i?ig2M"^9lp ~ (91 ' 92)^131/^1/^^ = , (36) 
(91 • q2)£q2tJ.iyp — £qiq2iypQltJ. ~ £qiq2tJ.pl2v — . (37) 

The two conditions reduce the set of independent covariants to two. This is the appropriate number of independent 
covariants since the photon transition is described by two independent amplitudes as e.g. by the El and AI2 transition 
amplitudes. 

Using the two constraint Eqs. (p6t and (|37|) the expansion (|34|) can be written in the form 



T^p. = (wi + W3- -^^^w}j e,i,2MP92. + (wi + M^2 - + (^J^^ £q,q2.pqip- (38) 

By comparing with the corresponding expressions in the Appendix one notes that the first and second terms in (j38p 
describe transitions into the longitudinal and transverse components of the J/ip. 
The quantities Wi are represented by the four-fold integrals 

oo 1 

W, = JdtJd''(3F,{t,l3i,P2,p3), (39) 



where we have suppressed the additional dependence of the integrand Fi on the set of variables p^, q^, ql; Wq, TOc, sx, sj/^ 
with sx = 1/A|^ and sj/^ = The integrals in Eq. ((39)) have branch points at — 4:{mq + rric)^ [diagram 

in Fig. [T]-a] and at = Arn^ [diagrams in Figs. [l]-b,c,d]. At these points the integrals become nonanalytical in the 
conventional sense when t — > oo. In order to check on the gauge invariance of the amplitude T^p^iQi, Q2)j we have 
taken the X-meson momentum squared to be below the closest unitarity threshold, i.e. < 4m^. We have checked 
explicitly that, for mx = 3.1 GeV and Wj/^ = 2.9 GeV, the gauge condition W4 = —W5 is numerically satisfied to 
very high accuracy. Note that the gauge invariance condition is independent of the overall couplings gx and gj/^ and 
thus the numerical check can be done irrelevant of their values. 

In the next step we introduce an infrared cutoff on the upper limit of the t-integration in Eq. (j39p . In this 
manner one removes all possible nonanalytic structures and thereby one obtains entire functions for the amplitudes, 
i.e. one has effectively instituted quark confinement, see Refs. mm. The value of A = 181 MeV was found by fitting 
the calculated basic quantities to the experimental data. However, for such a value of A the contributions coming 
from the bubble diagrams in Figs. [T]-b,c,d blow up at p^ = m\ compared with the contribution from the diagram 
Fig. [T]-a. The bubble diagrams are needed only to guarantee the gauge invariance of the matrix element. For physical 
applications one should take into account only the gauge invariant part of the diagram Fig. [T]-a. 

It is convenient to present the decay width via helicity or multipole amplitudes. The projection of the Lorentz 
amplitudes to the helicity amplitudes is given in Appendix. One has 

r,X^..«.^i|(|H.|» + |H.|^) . J-i|(M..PH-M,„P). (40) 



where the helicity amplitudes Hl and Ht are expressed in terms of the Lorentz amplitudes as 



W1+W3 ^W^4 

"^x|92| 



Ht = -imx\q2\ 



W1 + W2- (1 + — ) W4 



mx\q2\ 



9 9 



|92| = ^V^- (4^) 
2mx 
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The El and M2 multipole amplitudes are obtained via Aei/m2 = {Hl T Ht)/V^ ■ If we choose Ax = 3.0 GeV for 
the size parameter of the X(3872) we obtain Am2/Aei = 0.11, i.e. the electric multipole amplitude Aei dominates 
the transition, as expected. Nevertheless our predicted angular decay distribution W{{)) ^ 1 — 0.52cos^i? differs 
noticeably from its form W{'d) ^ 1 — 0.333 cos^i? for El dominance. It would be interesting to experimentally check 
on this prediction of our tetraquark model. 

In Fig. [2] we show a plot of the size parameter dependence of the decay width r{Xi — > J/ip + 7) together with the 
decay width r{Xi — ?> J/ip + 2n) taken from [l|. We correct an error of Ref. [ij in the normalization condition of the 
X meson, which led to a < 30% underestimate of the strong decay widths. Both decay widths become smaller as 
the size parameter increases. Note that the radiative decay width for Xh ~ —XuSinO + XdCosO is almost an order 
of magnitude smaller than that for Xi = XuCOsO + Xds'mO. If one takes Ax G (3,4) GeV with the central value 
Ax =3.5 GeV then our prediction for the ratio of widths reads 



T{Xi ^ 7 + J/V-) 



r{Xi J/ij + 27r 

which fits very well the experimental data from the Belle Collaboration [l 



= 0.15 ±0.03 (42) 

thcor 



0.14 ±0.05 Belle [10] 

(43) 



T{X J/ij 2tt) [ 0.22 ± 0.06 BABAR [13] 
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FIG. 2. The dependence of the decay widths T{Xi — 7 + J/tp) and r{Xi — >■ J/tp 2tt) on the size parameter Ax- 



V. SUMMARY AND CONCLUSION 



We have used our relativistic constituent quark model which includes infrared confinement in an effective way 
to calculate the radiative decay A" — 7 + J/tp. We take the X(3872) meson to be a tetraquark state with the 
quantum numbers J^'~" — l'^'^ . In order to introduce electromagnetic interactions we have gauged a nonlocal effective 
Lagrangian which describes the interaction of the X(3872) meson with its four constituent quarks by using the P- 
exponential path-independent formalism. We have calculated the matrix element of the transition A" — ^ 7 + J/tp and 
have shown its gauge invariance. We have evaluated the A' ^ 7 + J/ip decay width and the polarization of the J/ip in 
the decay. The calculated decay width is consistent with the available experimental data for reasonable values of the 
size parameter of the X(3872) meson. 
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Appendix A: Helicity and multipole amplitudes 

The material presented in this Appendix is adapted from similar material written down in in a slightly different 
context. There are two independent helicity amplitudes H\^-\^\j^^ which we denote by Hi (i = L,T) according to 
the helicity of the final meson state J/ip, where Xj/^ = and Xj/^ — ±1 stand for the longitudinal and transverse 
helicities of the J/tp. From parity one has if+;-o = ^-ff-;+o = Hl and -ffo;++ = — -^O; = Ht- 

We seek a covariant representation for the longitudinal and transverse projectors P^''^ which, when applied to the 
transition amplitude 7)ip;/, project onto the helicity amplitudes Hl^t according to 

iJ, ^Pfp^^,, {i = L,T). (Al) 

The projectors are defined by 

ipr = I (4(+)4'(-) - 4(-)£7^^(+)) 4/;(o) ' 

= ^ 4(0) {el''{+)slj,{+) - e-YHelj^i-)) , (A2) 

where we use the Jacob- Wick convention for the helicity polarization four-vectors as written down in [25|. The 
^-direction is defined by the momentum of the J/ip. The bars in the polarization four-vectors ef!^{Xj) of the photon 
are a reminder that the photon helicities are defined relative to the negative z-direction. In the present context it is 
important to take into account both parity configurations related by a helicity reflection in the definition of Eq. (|A2I) . 
In explicit form one has in the X-rest frame 



ex^(±) = (0;±l,j,0), K = (mx;0,0,0), 
exM(O) = (0;0,0,-l), 

4/,.(0) = ;;i^(l'r2l;0,0,-4±^), 

£-tp(±) = ^(0;Tl,-*,0), <Z2" = |g2|(l;0,0,-l). (A3) 
A convenient covariant representation of the projectors can be obtained in the form 

where 



hr - ^^^^ e'^"'^^^ q!^ , = -'oT^^ q^eP""''' , (A5) 

2 [qi ■ q2Y 2 (qi • 



and where the massive propagator functions are given by {V = X, J/ip) 



Sv a'iPv) = -9 a' + —-2 ■ (^6) 



The massive propagator functions are needed in the projectors Eq. (jA4p to project out the appropriate three- 
dimensional subspaces in the respective rest systems of the spin 1 particles. For the photon one exploits the gauge 
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freedom to write the propagator function as {—g^pi)- Note that the compact form (|A4|) is only obtained if one uses 

2U1J. 



the summed form (|A2p . The projection operators are orthonormal in the sense that P^'"'P|^ 
The angular decay distribution in the decay X{3872) — J' 7 + J/ip{-^ £^£~) is given by 



dV 



dcos i3 



= BR{J/^ ^ e+n^ || (J sin^ ^ + ^(1 + cos^ ^) \Ht\') , (A7) 



where d is the polar angle of either of the leptons £^ relative to the original flight direction of the J/tp, all in the rest 
system of the J/0. 

One can alternatively describe the transition in terms of the two multipole amplitudes Aei and Am2- The multipole 
amplitudes are related to the helicity amplitudes via i26i] 



V2 



M2 = 1 Hl + Ht 



(A8) 



The corresponding projectors onto the multipole amplitudes are given by 



"72 



^ A/2 



(A9) 



In Table 1 we have summarized the helicity and multipole amplitudes resulting from the relevant projections of 
the basic covariants Eq. ([5^ . The entries can be seen to satisfy the constraint equations Eqs. p6l37p . The multipole 
amplitudes Aei,m2 calculated from the gauge invariant structures Kj^pi, (i — 2,4,5,6) show the appropriate lowest- 
order power behavior Aei ~ I92I and Am2 ^ \Q2\'^- 

The leading \q2\ contribution to the angular decay distribution proportional to I^BiP is thus given by W{cos'd) cx 
(3 — cos^i?). The next-to-leading contribution proportional to 2TI{AeiA*j^j2) is down by one power of \q2\. The 
nonleading angular distribution is given by P^(cos I?) cx (1 - 3 cos^ ■&) (in the same units). 



4(0 _ Tpppf jy-y 
■"-M2 — M2 J^p.pi' 



Qlll£piyqiq2 
Qlp£pi^qiq2 
Q2l^£lJ.pqiq2 



kZI - {qiq2) Kill 








■ mv I -* I 



-'imxmj/^\q2\ 



tmx\q2\ 





-i\q2 



7M I ^ I 



V2 2mx 



^rnxm 



V2 mj/^ 
- mx +rn 



V2 

V2 2 



q2\ 



V2 {rnx+mj/^) 



^\q2? 



•J2 raji^ 
2m v 



mi 



V2 "ij/,p{mx+m,^,p) m\ 



2m\ 



V2 {rnx+mj/^,) 



^\q-A' 



TABLE I. Helicity and multipole projections of the basic tensors K^^pi^. The tensors K^l^, {i = 2,4,5,6) are gauge invariant. 
They satisfy g^if^p'^ = 0. 



[1] S. Dubnicka, A. Z. Dubnickova, M. A. Ivanov, J. G. Korner, "Quark model description of the tetraquark state X(3872) 
in a relativistic constituent quark model with infrared confinement", Phys. Rev. D81, 114007 (2010), arXiv: 1004. 1291 
[hep-ph]]; S. Dubnicka, A. Z. Dubnickova, M. A. Ivanov, J. G. Koerner, G. G. SaiduUaeva, AIP Conf. Proc. 1343, 385-387 
(2011). (arXiv:1 011.4417 [hep-ph]]. 

[2] T. Bran z, A. Faessler, T. Gutsche, M. A. Ivanov, ,1. G. Korner and V. E. Lyubovitskij, Phys. Rev. D 81, 034010 (2010), 
;arXiv:0912.3710 [hep-ph]]. 

[3] L. Maiani, F. Piccinini, A. D. Polosa and V. Riquer, Phys. Rev. D 71, 014028 (2005), jarXiv:hep-ph/0412098 . 
[4] S. K. Choi et al. [Belle CoUaboration], Phys. Rev. Lett. 91, 262001 (2003). 
[5] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 



13 



B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 93, 041801 (2004). 

D. E. Acosta et al. [CDF Collaboration], Phys. Rev. Lett. 93, 072001 (2004). 

V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 93, 162002 (2004). 

T. Aaltonen, et al. [CDF Collaboration], Phys. Rev. Lett. 103, 152001 (2009) a rXiv: 0906 .5218 ': Kai Yi, "Heavy quark 

meson spectroscopy at CDF", arXiv:0910.3163 . 

K. Abe et al. [ Belle CoUaboration ], hep-ex/0505037 . 

B. Aubert et al. [ BABAR Collaboration ], Phys. Rev. D74, 071101 (2006). hep-ex/0607050] . 
G. Gokhroo et al. [Belle CoUaboration], Phys. Rev. Lett. 97, 162002 (2006). 

E. Klempt and A. Zaitsev, Phys. Rept. 454, 1 (2007) arXiY;0708,.4 016 [hep-ph]]. 

B. Aubert et al. [ BABAR Collaboration ], Phys. Rev.' Lett. 102, 1 32001 (2009) arXiv:0809.0042l [hep-ex]]. 
T. Barnes and S. Godfrey, Phys. Rev. D 69, 054008 (2004) arXiv : hep- p h/03 1 1 1 62] . 

E. S. Swanson, Phys. Lett. B 598, 197 (2004) arXiv:hep-ph/0406080 . 

Y. -b. Dong, A. Faessler, T. Gutsche et al., Phys. Rev. D77, 094013 (2008) arXiv:080 23610l [hep-ph]]. 
M. Nielsen, C. M. Zanetti, Phys. Rev. D82, 116002 (2010) arXiv: 1006.0467 [hep-ph]]. 
M. A. Ivanov, M. P. Locher and V. E. Lyubovitskij , Few Body Syst. 21, 131 (1996). 
S. Mandelstam, Annals Phys. 19, 1 (1962); J. Terning, Phys. Rev. D 44, 887 (1991). 

A. Salam, Nuovo Cim. 25, 224 (1962); S. Weinberg, Phys. Rev. 130, 776 (1963); for review, see K. Hayashi et al.. Fort. 
Phys. 15, 625 (1967). 

G. V. Efimov and M. A. Ivanov, The Quark Confinement Model of Hadrons, (lOP Publishing, Bristol & Philadelphia, 
1993). 

J. G. Korner, M. C. Mauser, Lect. Notes Phys. 647, 212-244 (2004). 'hep-ph/0306082l. 
J. G. Korner, J. H. Kiihn, M. Krammer and H. Schneider, Nucl. Phys. B 229 (1983) 115. 
P.R. Anvil and J.J. Brehm, Phys. Rev. 145 (1966) 1152. 
W. N. Cottingham and B. R. Pollard, Annals Phys. 105 (1977) 111. 



